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ABSTRACT 

The '"OH/'^'OH and OD/OH isotope ratios are measured in the Oort-Cloud comet C/2002 T7 (LINEAR) through ground-based 
observations of the OH A^E+ - X^Ilj ultraviolet bands at 3063 A (0,0) and 3121 A (1,1) secured with the Very Large Telescope 
(VLT) feeding the Ultraviolet- Visual Echelle Spectrograph (UVES). From the '*OH/"*OH ratio, we find '«0/"*0 = 425 ± 55, equal 
within the uncertainties to the terrestrial value and to the ratio measured in other comets, although marginally smaller. We also estimate 
OD/OH from which we derive D/H = 2.5 ± 0.7 lO""* in water. This value is compatible with the water D/H ratios evaluated in other 
comets and marginally higher than the terrestrial value. 
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1. Introduction 

The determination of the abundance ratios of the stable isotopes 
of the light elements in different objects of the Solar System pro- 
vides important clues into the study of their origin and history. 
This is especially true for comets which carry the most valuable 
information regarding the material in the primitive solar nebula. 

The '^0/'**0 isotopic ratio has been measured from space 
missions in a few comets. In-situ measurements with the neu- 
tral and ion mass spectrometers onboard the Giotto spacecraft 
gave 16Q/18Q = 495+37 for H2O in comet IP/Halley (Eberhardt 
et al. 1199 5V A deep integration of the spectrum of the bright 
comet 153P/2002 CI (Ikeya-Zhang) with the sub-millimeter 
satelUte Odin led to the detection of the H^^O Une at 548 GHz 
(Lecacheux et al. 2003). Subsequent observations resulted in 
the determination of i''0/'**0 = 530+60, 530+60, 550+75 and 
508+33 in the Oort-Cloud comets Ikeya-Zhang, C/20 01 Q4 , 
C/2002 T7 and C/2004 Q2 respectively (Biver et al. l 20D71 l. 
Within the eiTorbars, these measurements are consistent with the 
terrestrial value ("' 0/"^0 ( SMOWQ) = 499), although marginally 
higher (Biver et al. I2007I I. More recently, laboratory analyses 
of the silicate and oxide mineral grains from the Jupiter family 
comet 8 IP/Wild 2 returned by the Stardust space mission pro- 
vided '*0/''^0 ratios also in excellent agreement with the teiTes- 
trial value. Only one refractory grain appeared marginally de- 
pleted in ''^O ("'0/"*0 = 576+78) as o bserve d in refractory in- 
clusions in meteorites (McKeegan et al. l20061 l. 

The D/H ratio has been measured in four comets. In-situ 
measurements provided D/H = 3.16+0.34 lO"'* for H2O in 
IP/Halley (Eberhardt et al. T99T. Balsiger et al. T^97), a fac- 
tor of two higher than the terrestrial value (D/H (SMOW) = 1 .556 
10 The advent of powerful sub-millimeter telescopes, namely 
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the Caltech Submillimeter Observatory and the James Clerck 
Maxwell telescope located in Hawaii, allowed the determina- 
tion of the D/H ratio for two exceptionally bright comets. In 
comet C/1996 B2 (Hyakutake), D/H wa s foun d equal to 2.9+1.0 
10 '* in H2O (Bockelee-Morvan et al. |1998l l, while, in comet 
C/1995 01 (Hale-Bopp), the ratios D/H = 3.3+0.8 lO'^^ in H2O 
and D/H = 2.3+0.4 10"^ in HCN were measured (Meier et al. 
1998a. 1998b ), confirming the high D/H value in comets. Both 
Hyakutake and Hale-Bopp are Oort-Cloud comets. Finally, bulk 
fragments of 8 IP/Wild 2 grains returned by Stardust indicated 
moderate D/H enhancements with respect to the terrestrial value. 
Although D/H in 8 IP/Wild 2 cannot be ascribed to water, the 
measured values overlap the range of water D/H ratios deter- 
mined in the other comets (McKeegan et aI. I2006l) . 

Among a series of spectra obtained with UVES at the VLT 
to measure the '"^N/'^N and '^C/'^'C isotope ratios in various 
comets from the 3880 A CN ultraviolet band (e.g. Arpigny et 
al. |2 003l H utsemekers et al.'2005^ Jehin et al. 2006 Manfroid 
et al. I2OO8 I1, we found that the spectrum of C/2002 T7 appeared 
bright enough to detect the "*OH Hnes in the A ^1+ -X ^11, bands 
at 3100 A allowing -for the first time- the determination of the 
i6q^isq j-^fjQ from ground-based observations. We also realized 
that the signal-to-noise ratio of our data was sufficient to allow a 
reasonable estimate of the OD/OH ratio from the same bands. 

The possibility of determining the '*0/'^0 ratio from the 
OH ultraviolet bands has been emphasized by Kim (20001 1. 
Measurements of the OD/OH ratio were already attempted by 
A'Hearn et al. ( |1985l l using high resolution spectra from the 
International Ultraviolet Explorer and resulting in the upper limit 
D/H < 4 10""* for comet C/1989 CI (Austin). These observations 
now become feasible from the ground thanks to the high ultravi- 
olet throughput of spectrographs like UVES at the VLT. 



2. Observations and data analysis 

Observations of comet C/2002 T7 were carried out with UVES 
mounted on the 8.2m UT2 telescope of the European Southern 
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Table 1. Observing circumstances 



Date 


r 


r 


A 


Offset 




Airmass 


(2004) 


(AU) 


(km/s) 


(AU) 


(10^ km) 


(s) 




May 6 


0.68 


15.8 


0.61 


1.3 


1080 


2.2-1.9 


May 26 


0.94 


25.6 


0.41 


0.0 


2677 


1.3-1.8 


May 26 


0.94 


25.6 


0.41 


0.0 


1800 


2.1-2.7 


May 28 


0.97 


25.9 


0.48 


10.0 


3600 


1.3-1.7 



r and r are thie comet lieliocentric distance and radial velocity; A is 
the geocentric distance; t is the exposure time; Airmass is given at the 
beginning and at the end of the exposure 



Observatory VLT. Spectra in the wavelength range 3040 A- 
10420 A were secured in service mode during the period May 6, 
2004 to June 12, 2004. The UVES settings 346-1-580 and 
437-H860 were used with dichroic #1 and #2 respectively. In 
the following, only the brighest ultraviolet spectra obtained on 
May 6, May 26 and May 28 are considered. The 0.44 x 10.0 
arcsec slit provided a resolving power R ^ 80000. The slit was 
oriented along the tail, centered on the nucleus on May 26, and 
off-set from the nucleus for the May 6 and May 28 observations. 
The observing circumstances are summarized in Table [T] 

The spectra were reduced using the UVES pipeline 
(Ballester et al. I2000I I, modified to accurately merge the orders 
taking into account the two-dimensional nature of the spectra. 
The flat-fields were obtained with the deuterium lamp which is 
more powerful in the ultraviolet. 

The data analysis and the isotopic ratio measurements were 
performed using the method designed to estimate the carbon 
and nitrogen isotopic ratios from the CN ultraviolet spectrum 
(Arpigny et al. ,2003. Jehin et al. .2004 and Manfroid et al. 
I2005I I. Basically, we compute synthetic fluorescence spectra of 
the i^OH, i^OH and '^^OD for the A - X^n, (0,0) and (1,1) 
ultraviolet bands for each observing circumstance. Isotope ratios 
are then estimated by fitting the observed OH spectra with a lin- 
ear combination of the synthetic spectra of the two species of 
interest. 

2A. The OH model 

We have developed a fluorescence model for OH similar to the 
one described by Schleicher and A'Hearn ( I1988I I. As lines of 
the OH(2-2) bands are clearly visible in our spectra we have 
included vibrational states up to v = 2 in the A^Z^ and X^fl, 
electronic states. For each vibrational state rotational levels up 
to 7 = 11/2 were included, leading to more than 900 electronic 
and vibration-rotation transitions. The system was then solved 
as described in Zucconi and Festou ( 119851 ). 

Accurate OH wavelengths were computed using the spectro- 
scopic constants of Colin et al. (2002) and Stark et al. (1994). 
OD wavelengths were computed using the spectroscopic con- 
stants of Abrams et al. (1994) and Stark et al. (1994). '^OH 
wavelengths were derived from the '^OH ones using the standard 
isotopic shift formula; they are consistent with the measured val- 
ues of Cheung et al. ( |19951 l. 

Electronic transition probabilities for OH and OD are given 
by Luque and Crosley ( 119981 11999I I. We used the dipole mo- 
ments of OH and OD measured by Peterson et al. (|1984| l to com- 
pute the rotational transition probabilities and the vibrational 
lifetimes computed by Mies (1974). Because of the very small 
difference in the structure of '^OH and '^OH the transition prob- 
abilities for "*0H and '^OH are the same. 



The OH fluorescence spectrum is strongly affected by the 
solar Fraunhofer lines, especially in the 0-0 band, so a care- 
fully calibrated solar atlas is required. We have used the Kurucz 
( |2055l ) atlas above 2990 A and the A'Hearn et al. ( lT983] l atlas 
below. 

The role of collisions in the OH emission, in particular those 
with charged particles inducing transitions in the A doublet 
ground rotational state, was first pointed out by Despois et al. 
( 119811 ) in the context of the 18 cm radio emission and then also 
considered in the UV emission by Schleicher ( 1983 Schleicher 
and A'Hearn 1988). Modeling the effect of collisions may be 
done by adding the collision probability transition rate between 
any two levels, / and j: 

= ^ n,{Y) v,.(r) cr,(/, v^.) 

c 

where the sum extends over all colliders, is the local density of 
the particles inducing the transition, v^. is the relative velocity of 
the particles and cr^. is the collision cross section. It also depends 
on the energy of the collision i.e. of v^.. The reciprocal transition 
rates are obtained through detailed balance: 

Cj,i = Cij—exp(Eij/kT) 
8j 

in which gi is the statistical weight and Ejj is the energy separa- 
tion between the states. In order to reduce the number of param- 
eters required to model the collisions we have adopted a simpli- 
fied expression of the form Cij = q\ for the transition in the 
A doublet ground state. In order to better fit the OH spectra we 
have also found necessary to take into account rotational excita- 
tion through a similar expression C/j = qni with different 
from only for dipole transitions, i.e. when AJ < 2, which ap- 
peared to correctly fit the data. Furthermore, since OH and OD 
have similar dipole moments, we assumed that collisional cross- 
sections are identical for both molecules. 

The model assumes that the '^OH lines are optically thin. 
This is verified by the fact that it correctly reproduces both the 
faint and strong OH emission lines. 

2.2. ^^OHA^OH 

Two i**OH lines at 3086.272 A and 3091.046 A are clearly 
detected in the (0,0) band. However these lines are strongly 
blended with the ~ 500 times brighter '^OH emission lines and 
then not useful for an accurate flux estimate. In fact the (1,1) 
band at 3121 A, while fainter, is better suited for the determina- 
tion of '''OH/'**OH since (i) the wavelength separation between 
"*OH and ""OH is larger (^ 0.3 A instead of 0.1 A), and (ii) 
the sensitivity of UVES rapidly increases towards longer wave- 
lengths while the atmospheric extinction decreases, resulting in 
a better signal-to-noise ratio. 

Fig. [1] illustrates a part of the observed OH (1,1) band to- 
gether with the synthetic spectrum from the model. Two "^OH 
lines are clearly identified. 

To actually evaluate '^OH/'^OH we first select the 3 brighest 
and best separated ''^OH lines at A ^ 3134.315 A, 3137.459 A 
and 3 142.203 A. These lines are then doppler- shifted and co- 
added with proper weights to produce an average profile which 
is compared to the ""OH profile similarly treated (cf. Jehin et al. 
12004 for more details on the method). We verified that the '^OH 
faint wings and nearby prompt emission lines (analysed in de- 
tail in a forthcoming paper) do not contaminate the '^'OH lines 
nor the measurement of the isotopic ratios. The ratio "'OH/"*OH 
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Fig. 1. A section of the May 6 spectrum of the OH (1,1) band 
in comet C/2002 T7. Black: observed spectrum; blue: synthetic 
emission spectrum of ^^OH fitted to the observed one (the ex- 
tended wings of the '^OH lines are not modelled here; on the 
other hand, in this and all three other figures, the presence of 
prompt emission is taken into account); red: synthetic fluores- 
cence spectrum of '^OH with the typical ratio '''OH/'^^OH - 500. 
The position of the "^OH lines at 3134.315 A and 3137.459 A is 
indicated. 
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Fig.2. Co-addition of the i**OH lines at ^ = 3134.315 A, 
3137.459 A and 3142.203 A from the May 6 spectrum after 
proper wavelength shifts. Co-added '^OH is shifted to A=0 while 
^^OH appears at 0.3. Black: observed spectrum; red: syn- 
thetic combined spectrum of '^OH, '^OH with ^^OH/"*OH = 
410, and an empirical fit to the OH extended wings. 



is then derived through an iterative procedure which is repeated 
for each spectrum independently. For the spectra of May 6, 26 
and 28 we respectively derive i*'OH/"^OH = 410±60, 510±130 
and 380+290. The uncertainties are estimated from the co-added 
spectra by considering the rms noise in spectral regions adjacent 
to the ^^OH lines, and by evaluating errors in the positioning of 
the underlying pseudo-continuum (i.e. the dust continuum plus 
the faint wings of the strong lines). The weighted average of all 
measurements gives "'OH/'^^OH = 425 ± 55. 

Since OH is essentially produced from the dissociation of 
H2O, ""OH/'^OH represents the ""O/'^O ratio in cometaiy wa- 
ter, with the reasonable assumption that photodissociation cross- 
sections are identical for Hi'^O and H2'^0. 



2.3. OD/OH 

The detection of OD lines is much more challenging since one 
may expect the OD lines to be a few thousand times fainter than 
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Fig. 3. Co-addition of 27 OD lines from the May 6 spectrum 
after proper wavelength shifts. Co-added OD is shifted to A-0. 
Black: observed spectrum; red: synthetic combined spectrum of 
OH, OD with OD/OH = 4 10 and an empirical fit to the OH 
extended wings. 




-0.4 



-0.2 



0.2 
Wavelength (A) 



0.4 



Fig. 4. Same as Fig.[3]for the May 26 spectrum. 



the OH lines. Fortunately, the wavelength separation between 
OD and OH (> 10 A) is much larger than between '^OH and 
'*'OH such that both the (0,0) and (1,1) bands can be used with 
no OD/OH blending (apart from chance coincidences). Since no 
individual OD lines could be detected, we consider the 30 brigh- 
est OD lines (as predicted by the model) for co-addition. After 
removing 3 of them, blended with other emission lines, an aver- 
age profile is built with careful Doppler-shifting and weighting 
as done for '^OH. Only our best spectra obtained on May 6 and 
May 26 are considered, noting that the (0,0) band -which domi- 
nates the co-addition- is best exposed on May 26 while the (1,1) 
band is best exposed on May 6, due to the difference in airmass. 
The resulting OD line profiles are illustrated in Fig. [3] and |4] and 
compared to a synthetic spectrum computed with OD/OH - 4 
10""*. OD is detected as a faint emission feature which is present 
at both epochs. From the measurement of the line intensities, we 
derive OD/OH = 3.3+1.1 lO '* and 4.1±2.0 10 for the spec- 
tra obtained on May 6 and 26 respectively. The weighted aver- 
age is OD/OH = 3.5±1.0 10 The difference in the lifetime of 
OD and OH (van Dishoeck and Dalgarno 1984) does not sig- 
nificantly affect our results since the part of the coma sampled 
by the UVES slit is two orders of magnitude smaller than the 
typical OH scale-length. The uncertainties on OD/OH were es- 
timated as for '^OH/'^OH. Possible errors on the isotopic ratios 
related to uncertainties on the collision coefficients were esti- 
mated via simulations and found to be negligible. Even in the 
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hypothetical case that colhsions differently affect OD and OH, 
errors are much smaller than the other uncertainties, as expected 
since the contribution of collisions is small with respect to the 
contribution due to pure fluorescence. 

To estimate the cometary D/H ratio in water, HDO/H2O must 
be evaluated. While the cross-section for photodissociation of 
HDO is similar to that of H2O, the production of OD+H is 
favoured over OH+D by a factor around 2.5 (Zhang and Imre 
119881 Engel and Schinke |T988] l. Assuming that the total branch- 
ing ratio for HDO -> OD H- H plus HDO -> OH -H D is equal to 
that of H2O ^ OH -I- H, we find HDO/H2O ^ 1 .4 OD/OH. With 
D/H = 0.5 HDO/H2O, we finally derive D/H = 2.5+0.7 10 
in cometary water. The factor (OD-i-H)/(OH-i-D) - 2.5 adopted 
in computing the branching ratios for the photodissociation of 
HDO is an average value over the spectral region where the 
cross-sections peak. In fact (OD-i-H)/(OH-i-D) depends on the 
wavelength and roughly ranges between 2 and 3 over the spectral 
regions where absorption is significant (Engel and Schinke ri9881 
Zhang et al. 1 19891 Yi et al. 120071 ). Fortunately, even if we adopt 
the extreme ratios (OD-hH)/(OH-hD) = 2 or (OD+H)/(OH-i-D) = 
3 instead of 2.5, the value of the D/H isotopic ratio is not changed 
by more than 6%. 

3. Discussion 

We have measured the oxygen isotopic ratio "'0/"^0 - 425 ± 55 
from the OH A - X^H, ultraviolet bands in comet C/2002 
T7. Although marginally smaller, our value do agree within the 
uncertainties with "'0/"*0 = 550 + 75 estimated from observa- 
tions by the Odin satelHte (Biver et al. l2007l l. with the "'O/'^O 
ratios determined in other comets, and with the terrestrial value 
(Sect.[B. 

To explain the so-called "oxygen anomaly" i.e. the fact that 
oxygen isotope variations in meteorites cannot be explained by 
mass-dependent fractionation, models of the pre-solar nebula 
based on CO self-shielding were proposed, predicting enrich- 
ments, with respect to the SMOW value, of '^O in cometary wa- 
ter up to '^0/'**0 -415 (Yurimoto & Kuramoto'2004). Recently, 
Sakamoto et al. (2007) found evidence for such an enrichment 
in a primitive carbonaceous chondrite, supporting self-shielding 
models. The value of ""O/'^O we found in C/2002 T7 is also 
marginally smaller than the terrestrial value and compatible 
with these predictions. On the other hand, the measurement of 
i^O/ii^O = 440 + 6 in the solar photosphere (Ayres et al. |2006l 
cf. Wiens et al. 120041 for a review of other, less accurate, mea- 
surements) indicates that solar ratios may deviate from the ter- 
restrial ratios by much larger factors than anticipated, requiring 
some revision of the models. More observations are then criti- 
cally needed to get an accurate value of '^0/"^0 in comets, as- 
suming that cometary water is pristine enough and can be char- 
acterized by a small set of representative values. Namely, if self- 
shielding is important in the formation of the solar system, it 
is not excluded that significant variations can be observed be- 
tween comets formed at different locations in the solar system, 
like Oort cloud and Jupiter-family comets. 

We also detected OD and estimated D/H = 2.5 ± 0.7 10 
in water. Our measurement is compatible with other values of 
D/H in cometary water and marginally higher than the terres- 
trial value (Sect.lT). Our observations were not optimized for the 
measurement of OD/OH (neither for '^OH/'^OH) and one of our 
best spectra was obtained at airmass ~ 2 with less than 20 min of 
exposure time for a comet of heliocentric magnitude m, ^ 5 (for 
comparison, comet Hale-Bopp reached nir ^ -1). All these ob- 
serving circumstances can be improved, including observations 



at negative heUocentric velocities to increase the OD/OH fluo- 
rescence efficiency ratio (cf. figure 1 of A'Hearn et al. 119851 ). 
This opens the possibility to routinely measure both the '^O/'^O 
and D/H ratios from the ground, together with the '^C/'^C and 
'"^N/'^N ratios, for a statistically significant sample of comets 
of different types (e.g. Oort-cloud, Halley-type, and hopefully 
Jupiter-family comets although the latter are usually fainter). 
The measurement of D/H is especially important since it allows 
to limit the contribution of comets to the terrestrial water, the 
high D abundance implying that no more than about 10 to 30% 
of Earth's water can be attributed to comets (e.g. Eberhardt et 
al. ri995[ Dauphas et al. 2000, Morbidefli et al. 2000). However, 
only a full census of D/H in comets could answer this question. 
In particular, if Jupiter-family comets, thought to have formed in 
farther and colder places in the Solar System, are characterized 
by an even higher D/H, closer to the ratio measured in the inter- 
stellar medium water, then the fraction of cometary H2O brought 
onto the Earth could be even smaller. 
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